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Abstract 



We present the signatures and prospects for the indirect detection of a Dirac right- 
handed neutrino dark matter candidate in neutrino telescopes, cosmic positron exper- 
iments and gamma-ray telescopes. An example of such a dark matter candidate can 
be found in extra-dimensional models. In some constructions, Kaluza-Klein states 
with the gauge quantum numbers of a right-handed neutrino can have sizable gauge 
interactions with Standard Model particles. For instance, in 5D warped Grand Unified 
Theories, it has been shown that a Kaluza-Klein right-handed neutrino may be stable 
and otherwise a phenomenologically viable dark matter candidate. We find that the 
prospects for the indirect detection of such a WIMP are encouraging, particularly for 
neutrino telescopes and cosmic positron experiments. 



1 Introduction 



Despite our ignorance of the nature of dark matter, numerous techniques have been developed 
in an effort to detect it [1]. Indeed, if the dark matter consists of Weakly Interacting Massive 
Particles (WIMPs), as is suggested by our theoretical prejudices, it will have a tiny but non- 
vanishing probability of interacting with the ordinary matter of the Standard Model (SM). 

In addition to collider searches, direct and indirect detection techniques have been devised 
to search for particle dark matter. In the case of direct detection, one attempts to measure 
the recoil energy of target nuclei in an underground detector due to its elastic scattering with 
a WIMP. Thousands of WIMPs are expected to cross each square centimeter every second 
at the surface of the Earth. We could potentially detect up to a few tens of events per year 
at near future detectors, depending on the precise interactions of the WIMP. 

The goal of indirect detection techniques on the other hand, is to detect a flux of cos- 
mic rays resulting from the annihilation of WIMPs in regions such as the Galactic center, 
the Galactic halo, or the interior of the Sun. These annihilations can produce potentially 
observable fluxes of gamma-rays, neutrinos, positrons, anti-protons and anti-deuterons. 

Given the abundance of experimental activity related to dark matter detection, it is 
timely to study the distinctive signatures expected in different dark matter scenarios. Such 
signatures and event rates can vary substantially from one WIMP model to another. The 
most studied candidate so far is the Lightest Supersymmetric Particle (LSP) in models of Im- 
parity conserving supersymmetry. Some alternatives to neutralinos and other LSPs do exist, 
however. Lately, the possibility of Kaluza-Klein dark matter has been explored. While the 
idea that dark matter could be made of Kaluza-Klein (KK) particles is very tempting, it 
turns out that this is not so easy to achieve. Indeed, in most extra-dimensional models, there 
are no stable KK states, all being able to decay to SM particles. An exception is the class of 
models in Universal Extra Dimensions (UED) [2] . In this case, all SM fields propagate in flat 
toroidal extra dimensions. Translation invariance along an extra dimension is only broken 
by the orbifold imposed to recover a chiral SM spectrum. Still, there is a remnant discrete 
symmetry called KK parity, (—1)™, where n is the KK number. This symmetry insures that 
interaction vertices cannot involve an odd number of odd-KK states and, therefore, a vertex 
with two SM particles (with n — 0) and one KK state (with n — 1) is forbidden. As a 
result, the Lightest KK Particle (LKP) with n — 1 cannot decay into SM particles and is 
stable. For ~TeV _1 sized extra dimensions, the LKP can act as a WIMP. Relic density [3,4], 
direct [5,6] and indirect detection [5,7-12] studies of this candidate have all been carried out 
in the last few years. Constraints on these models from radion cosmology have also been 
studied [13]. 

The second example of KK dark matter arises in the context of warped geometries and 
more specifically in the context of warped Grand Unified Theories (GUTs) [14,15]. To locate 
these constructions in the landscape of extra dimensional models, recall that the interest in 
the phenomenology of extra dimensions over the last few years has been motivated by the 
goal of understanding the weak scale. The only extra-dimensional geometry which really 
addresses the hierarchy problem is the Randall-Sundrum geometry. Particle physics model 
building in this framework has been flourishing and a favorite class of models has emerged: 
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that where all SM fields propagate in the bulk of AdSs, except for the Higgs (or alternative 
physics responsible for electroweak symmetry breaking) which is localized on the IR brane. 
In addition, the electroweak gauge group should be extended to SU(2) L x SU(2) R x U(l). 
Those models were embedded into a GUT in Refs. [14, 15] and it is in this context that a 
viable dark matter appears. 

In these models, a stable KK fermion can arise as a consequence of imposing proton 
stability in a way very reminiscent to R-parity stabilizing the LSP in supersymmetric models. 
The symmetry is called Z3 and is a linear combination of baryon number and SU(3) color. 
It actually exists in the SM but SM particles are not charged under it since only colored 
particles carry baryon number in the SM. In Refs. [14,15], and more generally in higher 
dimensional GUTs, baryon number can be assigned in such a way that there exists exotic 
KK states with the gauge quantum numbers of a lepton and which carry baryon number 
as well as KK quarks which carry non-standard baryon number. These particles carry a 
non-zero Z 3 charge. The lightest of these, called the Lightest Z 3 Particle (LZP), is stable 
since it cannot decay into SM particles. 

The only massive elementary Dirac fermion (with a mass in the 1 GeV - 1 TeV range) 
which could be a viable dark matter candidate is the neutrino. If such a neutrino had the 
same coupling to the Z as in the SM, however, it would be excluded by direct detection 
experiments. Its coupling to the Z, therefore, must be suppressed 1 . Thus, we are left with 
the possibility of a KK Right-Handed (RH) neutrino. In models where the electroweak gauge 
group is extended to SU(2) L x SU(2) R x U(l), the RH neutrino has gauge interactions in 
particular with the additional Z'. Nevertheless, its interactions with ordinary matter are 
feeble because they involve the additional gauge bosons which have a large mass (> 3 TeV). 
This opens the possibility of a weakly interacting Dirac RH neutrino. 

In this work, we present the prospects for the indirect detection of right-handed neutrino 
dark matter, based on the particular example of the LZP. In section |2l we will present the 
properties of this particle, how it gets its mass and how it interacts with ordinary matter. 
Section |3] deals with the predictions for the detection of neutrinos generated in annihilation 
of the LZP in the Sun. Signatures due to cosmic positrons from annihilation in the galactic 
halo are presented in section 0] and gamma-rays from annihilation in the Galactic center in 
section^ We summarize our results and make comparisons with LSP and LKP dark matter 
candidates in the following section. 

2 A model of Dirac right-handed neutrino dark matter 

As mentioned in the introduction, the LZP arises as a viable dark matter candidate in warped 
GUTs. An interesting feature of these models is that GUT states such as X, Y gauge bosons 
appear at the TeV scale (via the KK excitations). In SO (10), there are also the X' , Y' , X s , 
Z' , etc. that the LZP can couple to. 

The LZP is a Kaluza-Klein fermion, which is a four-component spinor and vector-like 
object. It gets its mass by marrying the left-handed with the right-handed chirality. As 

1 Note that in SUSY, such constraints are weaker because of the Majorana nature of the neutralino. 
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explained in great detail in Ref. [15], it can be naturally very light, much lighter than the 
KK scale of Randall-Sundrum models, namely Mkk ^ 3 TeV. This is because the RH 
chirality is localized near the TeV brane while the LH one is near the Planck brane. The 
overlap of wave functions is small, resulting in a small Dirac KK mass. Its lightness is related 
to the top quark's heaviness but not entirely fixed by it, so that we will consider LZPs in 
the mass range of approximately 20 GeV to a few TeV. In this paper, we refer to the LZP 
as if it were a chiral fermion because only the RH chirality has significant interactions and 
the other chirality decouples. In addition, the LZP has the same gauge quantum numbers 
as the RH neutrino of SO (10) or Pati-Salam. As a result, we refer to it as a "Dirac RH 
neutrino" and denote it as u' R . 

The LZP couples to the TeV KK gauge bosons of 5*0(10). In addition, when electroweak 
symmetry is broken, Z — Z' mixing induces a coupling of the RH neutrino to the SM Z gauge 
boson. This coupling is suppressed by (M z / M z >) 2 ■ If M z > ~ few TeV (the mass of KK gauge 
bosons is set by Mkk), the size of this coupling will typically be ideal for a WIMP. There 
is actually a second source for the Z-LZP coupling, which we will not discuss here but refer 

v' 

the reader to Ref. [15] for a detailed explanation. We denote g z R as the coupling of v' R to 
the Z. There is enough freedom in the model under consideration to treat this coupling as 
an almost arbitrary parameter. 

As far as indirect detection is concerned, three phenomenological quantities are especially 
important: the annihilation cross section of the LZP, the fraction of annihilations which 
produce various SM particles and the LZP's elastic scattering cross section with nuclei (for 
the WIMP's capture in the Sun). These quantities were described in detail in Ref. [15]. 
For LZPs lighter than approximately 100 GeV, LZP annihilations proceed dominantly via 
s-channel Z-exchange. For larger masses, annihilation via the t-channel exchange of X s into 
top quarks or via s-channel Z' exchange into tt, bb, W + W~ and Zh dominates. 

Concerning elastic scattering, as is well-known for a Dirac neutrino, the Spin Independent 
(SI) elastic scattering cross section via a t-channel Z exchange has the form 



a si oc 



Z(l -4sin 2 fl w ) - (A-Z)} 2 



Since 4sm 2 6*w ~ 1, the coupling to protons is suppressed. In these models, the Spin- 
Dependent (SD) elastic scattering cross section can also play an important role, although 
only through scattering off of hydrogen in the Sun, since the spin of 4 He is zero. 

As presented in Ref. [15], there is a large region of parameter space in these models. The 

v' 

value of g z R depends sensitively on where we sit in this parameter space. In this study, we 
focus only on the regions in which the total annihilation cross section leads to a thermal relic 
density which does not exceed the measured value. It is indeed conceivable that the relic 
density we observe today is not only due to the contribution obtained from the standard 
thermal calculation, but also from non-thermal processes, such as the decay of heavier KK 
modes, gauge bosons or fermions. We have thus selected three representative regions of 
parameter space, corresponding to Mkk = 4, 6 and 10 TeV. For each Mkk value, we have 

v' 

two sets of values for g z R named 'minimal' and 'median'. The corresponding predictions for 
the dark matter relic density are given in Fig. and are summarized in Table ^ 
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nh2 M kk = 4TeV ^h 2 M kk = 6TeV &\\ 2 M kk = 10TeV 




Figure 1: Predictions for Q hZP h 2 as a function of the LZP mass for three values of the 
KK gauge boson mass, M KK . In each case, the upper dashed blue line corresponds to the 

v' 

'minimal' value of the Z-LZP coupling, g z R , while the lower pink curve is for a 'median' 
value of this coupling. In this article, we study the regions of parameter space given in the 
third column of Table 



3 Neutrinos From LZP Annihilations in the Sun 

WIMPs travelling through the Galactic halo will occasionally scatter off of large bodies, such 
as the Sun, and become gravitationally bound. The capture rate for LZPs in the Sun is given 
by [16]: 

n& AK -, a i8-i{ Piocai \ / 27Qkm/s \ 3 / Vh.sd + <r H ,si + 0-067 QHejA / 1000 GeV y 

6 " 4 - 5xlu 8 [oAGev/ctf) \^~r) [ io=^b 

(2) 

where pi OC ai and i>i OC ai are the local density and rms velocity of halo dark matter particles, 
respectively. <7h,sd, 0"h,si and o"H e ,si are the LZP's elastic Spin-Dependent (SD) and Spin- 
Independent (SI) scattering cross sections with hydrogen and helium nuclei. The cross 
section for LZP spin-independent scattering off of neutrons is much larger than for protons, 
making <7h,si unimportant compared to <7h,sd and o"H e ,si- 

The spin-independent cross sections in this model are [15]: 



9z 
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' e 2 /4(l - 4 sin 2 9 W ) 2 \?z R ) (2(1 - 4 sin 2 9 W ) - 2) 

aH ' SI = 64vrMf sin 2 6 W cos 2 6 W ' aHe ' SI = 64vrMf sin 2 6 W cos 2 6 W ' 

(3) 

where fJ,H,He = ^LZp^H,_H"e/ (w^lzp + rriH,He)- Note that aue,si dominates by a factor 10 4 over 

&H,SI- 

The cross section for the SD interaction is given by [15]: 

3 u 2 A 2 e q Vr 

cth,sd = , A = a p , a p = — 2 f , [-Au + Ad + As] , (4) 

7r 8M| cos # w sin W 
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Data sets 


Elastic scattering 
cross section, cr LZP _ nuc i eon 


LZP mass range 
leading to Vth 2 < 0.11 


M KK =4 TeV, Minimal gf 


~ 6 x 10~ 8 pb 


40 GeV<m LZP <45 GeV, 
m hZP > 700 GeV 


M KK =4 TeV, Median g% 


~ 2 x 10~ 7 pb 


35 GeV <m LZP < 50 GeV, 
m hZP > 350 GeV 


M KK =6 TeV, Minimal gfr 


~ 10~ 8 pb 


41 GeV < m LZP < 48 GeV, 
m LZP > 1500 GeV 


M KK =6 TeV, Median g v / 


~ 10~ 7 pb 


40 GeV<m LZP <50 GeV, 
m LZP >310 GeV 


M KK =10 TeV, Minimal gf 


~ 10~ 9 pb 


m LZP > 3100 GeV 


M KK =10 TeV, Median g£ 


~ 10~ 8 pb 


40 GeV, <m LZP <50 GeV 
m LZP > 900 GeV 



Table 1: Data sets used in our analysis (in particular in FigEJ). They can lead to the right 
relic density (see FigQJ) and do not exceed direct detection limits. 



where Au = 0.78 ±0.02, Ad = -0.48 ±0.02 and As = -0.15 ±0.02 [17] are the spins carried 
by the u, d and s quarks in the proton, respectively. 

The time rate of change of the number of LZPs in the Sun is given by: 

N = C G — A & N 2 , (5) 

where C & is the capture rate and A & is the annihilation cross section times the relative 
WIMP velocity per volume. A & is given by: 

A" = (6) 

where V e g is the effective volume of the core of the Sun and is determined by matching 
the core temperature with the gravitational potential energy of a single WIMP at the core 
radius. This was found to be [18]: 

2fi , / 1000 GeV\ 3/2 

V cS = 1.8 x 10 26 cm 3 . (7) 

V m LZP / 

The present day annihilation rate of LZPs in the Sun is given by: 

T = ^A Q N 2 = \c Q tanh 2 (VC°A^ t Q ) , (8) 
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where t Q ~ 4.5 billion years is the age of the solar system. If the capture rate in the Sun 
and annihilation cross section are large enough, these two processes will reach equilibrium 
with each other. This occurs for v C & A & t > 1. In fig. El V C & A Q t Q is evaluated at 
(av) k 3 x 10 -26 cm 3 /s, the value of the annihilation cross section leading to the correct 
relic density. Throughout the parameter region leading to the observed relic density, the 
Sun consistently reaches equilibrium between the LZP capture and annihilation rates. 




Figure 2: V C G A G t Q evaluated at (av) k 3 x 1CT 26 cm 3 /s, the value of the annihilation 
cross section leading to the correct thermal relic density. Throughout the parameter space, 
the Sun consistently reaches equilibrium between the LZP capture and annihilation rates. 
Figure from Ref. [15]. 

The next step is to determine the spectrum of neutrinos generated in the annihilation 
of LZPs. This depends on the region of parameter space considered. In low mass models 
(^lzp 100 GeV), LZPs annihilate primarily through s-channel Z exchange (via Z — Z' 
mixing), and therefore the annihilation modes very nearly follow the decay modes of the Z, 
ie. 6.7% to neutrino pairs of each flavor, 3.4% to charged lepton pairs of each flavor and the 
remaining ~70% to quark pairs. For more massive LZPs, annihilations to ti, W + W~ and 
Zh typically dominate (see figure 

For the first case (m^zp ^ 100 GeV), the neutrinos generated will be of sufficiently low 
energy that one could safely neglect their interactions with the solar medium. Since the 
experimental (muon) energy thresholds of high-energy neutrino telescopes are in the range 
of 10 to 100 GeV, few if any of the neutrinos produced in the decays of quarks or charged 
leptons produced will yield observable muons. 2 Focusing on annihilations directly to neutrino 

2 Neutrino telescopes detect muons generated in charged current interactions. These muons appear as 
tracks through the detector volume. Electrons produced by VeS simply generate an electromagnetic shower 
rather than an extended track. Showers can in principle be observed, but the energy thresholds for such 
events are a few TeV, so are not of interest to us here, r's decay too quickly to produce observable tracks. 
Only at PeV energies do r's travel further than ~ 100 meters. Again, decaying r's can produce showers, 
but well below the energy threshold. All neutrino flavors can produce hadronic showers via neutral current 
interactions. For our purposes here, only the z^v's and f e 's which have oscillated into v^s can be detected. 
Very few will come from v e because of the small corresponding mixing angle. 
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Figure 3: The fraction of LZP annihilations to various final states. Results are shown for 
friKK = 4 and 10 TeV in the left and right frames, respectively. 



pairs, the observed muon neutrino flux is a delta function at the mass of the LZP: 



dN Vlt 2TA uP S(m LZP ) 



dE Va 



(9) 



where T is the rate of LZP annihilations in the Sun, A vV m 



0.067 is the fraction of annihi- 
lations to Vp&y, (after including neutrino oscillations, A vD is the average of the annihilation 
fractions to muon and tau neutrinos) and -Res is the Earth-Sun distance. 

For heavier LZPs, annihilations to top quark pairs often dominate. Top quarks decay 
nearly 100% of the time to a W ± and b, each of which can generate neutrinos in their decays. 
The neutrino spectrum from this channel is given by [19]: 



dN u 



2YA f 



k dE v / tt 

where the first of these terms is given by 



(dN v 



\dE v 



t-*w ± 



(mf - mf^)7/9 



dN p 

Je }j 



In 



+ 



dN v 
dE„ 



(10) 



min(E i ,/( 7 (l-/?)),m t /2) 



for 



max(E^/( 7 (1 + /?)), m^/2m t ) 
7(1-/5) m 2 w /2m t < E v < 7 (1 + 0) m t /2, (11) 



where 7 = mizp/^t is the Lorentz boost of the produced top quark and (3 is the velocity 
corresponding to this Lorentz factor. „/± « 0.107 for each neutrino flavor. The second 
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term is given by: 

dN u \ B b ^ vl ± x m t 



dEJt^b if3Q.n{m 2 t -m 2 w ) 



3x 2 - 3y 2 + 4x 3 /3 - % 3 /3 + 2 In |) , 
for E v < 0.73 7 (1+ (3){m 2 -m 2 w )/2m t , (12) 



where x = 2m t E u /(0.73 7(1 + /3)(m t 2 + m^)) and y = min(l, 2m t E iy /(0.73 7(1 - [3)(m 2 + 
m^))). B b _, u i± x ~ 0.103 for electron and muon flavors and ~ 0.025 for tau flavor. 

After this spectrum of neutrinos is generated, it must be propagated through the Sun. 
Above E v ~ 100 GeV, this can be an important effect. The suppression of the neutrino flux 
as a result of absorption via charged current interactions in the Sun is given by [20] : 



dE v {Ev) ~ 


d$ ini 


x e 


GeV 

J 


d ^ (E ) ~ 
dE v [K) ~ 


iK <*> 


x e" 


-E„/200 GeV 


(E ) ~ 
dE v [Eu) ~ 




x e" 


-E„/130 GeV 


dE v [Ev) ~ 




x e 


-E v /200 GeV 



(13) 

where the spectra with the label, ini, are the initial neutrino spectra injected in the the core 
of the Sun. 

The effects of propagation on the tau neutrino spectrum are more complicated. Charged 
current interactions of tau neutrinos produce a tau lepton which, unlike muons or electrons, 
will decay generating another tau neutrino before coming nearly to rest. As a result of this 
regeneration process, the total number of tau neutrinos injected near the center of the Sun 
equals the number that escape, although many of these may have lost much of their energy. 
The unscattered component of the tau neutrino flux is similar to that of muon and electron 
neutrinos [20]: 



dE„ y ' dE, 



V 

ini 



w ^h* 7 ") X e"^ /230GeV (14) 

This is accompanied, however, by a spectrum of tau neutrinos which have been scattered 
in the solar medium. The distribution of the energies of the scattered component of tau 
neutrinos is roughly described by a log- normal distribution [20]: 



dN 1 

-n=r — ; — exp 

dE V^FlnlOaE 



2a 2 \Ej 



(15) 



with E t = 60 GeV and a = 0.53 for tau neutrinos and E t = 113 GeV and a = 0.49 for tau 
anti-neutrinos. 



8 



x) 1Q -5 



1 




l . . . lWUUU 



50 70 100 200 300 500 700 
E v (GeV) 



50 100 



200 500 1000 

E v (GeV) 



Figure 4: The spectrum of neutrinos and anti-neutrinos per LZP annihilation to top quark 
pairs in the Sun. The effects of solar propagation have been taken into account. Results 
in the left and right frames correspond to LZPs with masses of 800 and 1600 GeV, respec- 
tively. Light (green) dashed, light (green) dotted, dark (blue) dashed and dark (blue) dotted 
lines denote u^, v T and z/ r , respectively. The solid black line is the sum of these four 
contributions. We have used an annihilation cross section of 3 x!0~ 26 cm 3 /s. 



In figure |3] we show the final spectrum of neutrinos (plus anti- neutrinos) for LZP annihi- 
lations to ti. 

To detect a flux of high-energy neutrinos, large volume neutrino telescopes have been 
developed [21]. Muon neutrinos interacting in, or near, such an experiment generate energetic 
muons via charged current interactions which emit radiation as they propagate through the 
detector's Cerenkov medium, such as ice or water. The rate of neutrino-induced muons 
observed in a high-energy neutrino telescope is found by: 

^events " / / ^ R ^ ~ V) X E ») A ^ dE », ^ ( 16 ) 

where <t v n(E u ) is the neutrino-nucleon charged current interaction cross section [22], (1 — y) 
is the fraction of neutrino energy which goes into the muon, A e R is the effective area of the 
detector, -R M ((1 — y) x E v ) is the distance a muon of energy, (1 — y) E u , travels before falling 
below the muon energy threshold of the experiment, called the muon range [23]. 

Currently, the AMANDA telescope at the South Pole has been taking data for a number 
of years. With a muon energy threshold of approximately 25 GeV, AMANDA may be 
sensitive to LZPs in the lower 35 — 50 GeV. Using an effective area of 
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Figure 5: The number of events per year in a detector with a cubic kilometer effective 
area, such as IceCube. A 50 GeV muon energy threshold has been used. Solid (black) lines, 
dashed (blue) lines and dotted (red) lines correspond to Mkk =4, 6 and 10 TeV, respectively. 
For each value of Mkk, two lines are shown. These upper and lower lines use parameters 
corresponding to a 'median' and 'minimum' set of couplings (see Tabled for more details). 
Values are shown only over the mass ranges which do not overproduce the density of dark 
matter. 



50,000 square meters for the AMANDA experiment, we estimate an event rate of: 

N mna „ 2.9yr- ^». S D + ^. + 0.067^3. ^ 

for a LZP with a mass of 40-50 GeV. The ANTARES experiment, currently under construc- 
tion in the Mediterranean, is designed to have a somewhat lower energy threshold (~10 
GeV) and an 100,000 square meter effective area. For ANTARES, we estimate an event rate 

JV_ » 10.5 yr- ( gg±f5p + 0-067^ (lfj) 

for a LZP with a mass of 40-50 GeV. 

For heavier LZPs, which annihilate dominantly to tt, W + W~ or Zh, larger volume 
neutrino telescopes will most likely be needed. The IceCube experiment, currently under 
construction at the South Pole, will have a full cubic kilometer of instrumented volume. The 
event rate predicted in IceCube are shown in figure El 

These rates must compete with a background of atmospheric neutrinos [24] and neutrinos 
generated in cosmic ray interactions in the Sun's corona [25] . After cuts, approximately 5 to 
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10 atmospheric neutrino events are expected per year in AMANDA within the solid angle 
in the direction of the Sun. The background rate from cosmic ray interactions is expected 
to be significantly lower. The current data collected by AMANDA should be sufficient to 
exclude ~50 GeV LZPs with elastic scattering cross sections larger than ~ 10~ 6 pb, although 
a proper analysis of the data would need to be done to calculate a precise limit. 

The background rate of atmospheric neutrinos in ANTARES will be somewhat higher 
than the rate in AMANDA due to their larger area and lower energy threshold. When 
completed, ANTARES should be slightly more sensitive to ~50 GeV LZPs than AMANDA. 

IceCube, with a background rate of ~100 atmospheric neutrino events per year, will be 
able to test models which predict more than ~10 events per year (see figure EI). With several 
years of exposure, perhaps a rate of a few events per year could be detected. 



4 Signatures of LZP Annihilations in Cosmic Positrons 



High-energy positrons can be produced in LZP annihilations throughout the Galactic halo. 
These positrons then travel under the influence of tangled Galactic magnetic fields, losing en- 
ergy through inverse Compton scattering and synchrotron processes. The positron spectrum 
observed can be arrived at by solving the diffusion-loss equation: 



d dn e + 
dtdE P + 



V 



K(E e+ ,x)\7 



IE, 



+ 







dE P 



b{E e 



« dn e ^ 
,x)- 



dE, 



+ Q{E e+ ,x), 



(19) 



where dn e +/dE e + is the number density of positrons per unit energy, K(E e +,x) is the dif- 



fusion constant, b(E e +,x) is the rate of energy loss and Q{E e 
distribution of positrons injected through LZP annihilations. 
We parameterize the diffusion constant [26,27] by: 



x. 



is the spatial and energy 



K{E e+ ) 

and the energy loss rate by: 



3.3 x 10 



b(E e+ ) 



2S 



jO.47 



0.47 



cm 2 s , 



10" 16 £ 2 + s- 1 . 



(20) 



(21) 



b(E e +) is the result of inverse Compton scattering on both starlight and the cosmic microwave 
background [28]. We impose the boundary conditions for our diffusion zone corresponding 
to a slab of thickness 2L, where L is fit to observations to be approximately 4 kpc [26]. 
Beyond the boundaries of our diffusion zone, we drop the positron density to zero (free 
escape boundary conditions). These diffusion parameters are constrained by analysis of 
stable nuclei in cosmic rays (primarily by fitting the boron to carbon ratio) along with 
radioactive secondaries in the cosmic ray spectrum [29]. 

In addition to the propagation effects described above, as positrons approach the solar 
system, their interaction with the solar wind and magnetosphere can be important. These 
effects, called solar modulation, may considerably effect the positron spectrum at energies 
below 10-20 GeV. 
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The spectrum of positrons injected into the halo depends on the dominant annihilation 
modes of the LZP. For somewhat light LZPs (m^zp < 100 GeV), annihilations produce 
positrons directly (e + e~) as well as through the decays of muon, tau and quark pairs. The 
annihilations to e + e~, and to a lesser extent and t + t~, generate a very hard spectrum 

of positrons. Quark decays typically generate positrons with lower energies. For heavier 
LZPs, which annihilate to ti, W + W~ or Zh, a softer and less distinctive positron spectrum 
is generated. We have used PYTHIA [30], as implemented in DarkSusy program [31], to 
generate the injected positron spectra. 

Finally, the spatial distribution of positrons injected into the halo depends on the halo 
model assumed and on the quantity of dark matter substructures distributed in the local 
halo. Since positrons in the range of energies considered here typically lose most of their 
energy over distances greater than a few kiloparsecs, this channel does not efficiently sample 
annihilations near the Galactic center. For this reason, the presence of a cusp in the Galactic 
halo distribution is of only secondary importance. We have used a standard Navarro-Frenk- 
White halo profile in our calculations [32]. 

The presence of dark substructure in the local halo can lead to enhancements in the dark 
matter annihilation rate. Analytic estimates of this effect find that the positron flux is likely 
to be boosted by a factor of 2 to 5 due to such effects [33]. Recently, however, it has been 
suggested that greater quantities of substructure may have formed [34] and, depending on 
how many of these object survive [35], considerably larger boost factors may be possible. 

The HEAT experiment (High-Energy Antimatter Telescope) has made the most detailed 
measurements of the high-energy cosmic positron spectrum to date. This data appears 
to contain an excess in the spectrum, peaking at around 7-10 GeV and extending to higher 
energies [36] . Although the source of these positrons is not known, it has been often suggested 
that this signal could be the product of dark matter annihilations [9,27,37-39]. 

In figure El we show the positron spectrum from LZP annihilations compared to the 
HEAT data. The results are shown as the ratio of positrons to positrons plus electrons, 
called the positron fraction. Included are the primary and secondary electron fluxes as well 
as the secondary positron flux calculated in the model of Ref. [40]. The solid line in the 
figure is for this background alone, with no dark matter component. The dashed and dotted 
lines correspond to LZP masses of 50 and 40 GeV, respectively. These have been normalized 
using boost factors of 23 and 17, respectively, assuming a local dark matter density of p = 0.4 
GeV/cm 3 . These models both fit the HEAT data very well, yielding \ 2 ~ 1-1 P er degree of 
freedom. 

Although considerably heavier LZPs cannot produce the excess observed by HEAT with- 
out unnaturally large boost factors, these models may still be accessible to future cosmic 
positron experiments, which will be considerably more sensitive to dark matter annihila- 
tions [41]. 

The PAMELA experiment, to be launched into orbit later this year, is a satellite borne 
experiment designed to study the matter-antimatter asymmetry of the universe with very 
precise cosmic ray measurements. PAMELA'S primary objectives include the measurement 
of the cosmic positron spectrum up to 270 GeV [42], well beyond the range studied by HEAT. 
This improvement is made possible by PAMELA'S large acceptance (20.5 cm 2 sr [43]) and 
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Figure 6: The cosmic positron fraction from LZP annihilations for masses of 40 (dots) and 
50 (dashes) GeV. In each case, the normalization of the dark matter contribution was chosen 
to maximize the quality of the fit to the HEAT data (shown as error bars). We find x 2 's of 
approximately 1.1 per degree of freedom for each case. Using an annihilation cross section 
of 3 x 10~ 26 cm 3 /s, and a local dark matter density of p(local)=0.4 GeV/cm 3 , boost factors 
of 23 and 17 were needed to fit the data. The solid line is the background-only prediction. 



long exposure time (3 years). 

The AMS-02 experiment has been designed to be deployed on the International Space 
Station (ISS) for a three year mission sometime around the end of the decade. AMS-02's 
acceptance of about 450 cm 2 sr is considerably larger than PAMELA'S. This, along with su- 
perior energy resolution and electron and anti-proton rejection, makes AMS-02 the premiere 
experiment for measuring the cosmic positron spectrum. 

If LZP annihilations are in fact responsible for the HEAT excess, PAMELA and AMS-02 
should be able to unambiguously identify this fact. In the left frame of figure we show the 
spectrum of figure^ with the error bars projected for the PAMELA (larger, red) and AMS-02 
(smaller, blue) experiments. In addition to the overall precision of these measurements, the 
distinctive feature at the LZP mass would be unmistakable to these experiments. In the right 
frame of figure we show this same model, but with a conservative boost factor of 2. Even 
in this case, the distinctive spectral feature at the LZP mass can be confidently identified. 
It may also be possible to identify a discrete feature in the cosmic positron spectrum using 
Atmospheric Cerenkov Telescopes (ACTs) [11]. 

If LZPs are heavier and annihilate mostly to top quarks, or gauge and higgs bosons, 
positrons will be produced in the subsequent decay chains, i.e. t — > W ± b, W + — > e + z/, 
b — > e + uX, etc. The resulting spectrum, for a 400 GeV LZP annihilating to ti, is shown in 
figure |H1 Although no highly distinctive features are present in this spectrum, it may still 
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Figure 7: The cosmic positron fraction from 50 GeV LZP annihilations compared to the 
error bars projected for the PAMELA (larger, red) and AMS-02 (smaller, blue) experiments. 
In the left frame, the spectrum is normalized to the HEAT data (see figure ^ with a boost 
factor of 23. In the right frame, a conservative boost factor of 2 was used. The solid line is 
the background-only prediction. 



be possible to statistically identify such a spectrum over the background prediction. 

Other anti- matter channels for indirect dark matter searches have been studied as well, 
namely cosmic anti-proton and anti-deuteron studies. Although such searches are indeed 
interesting, they sample over considerably larger regions of the Galactic halo, making the 
choice of diffusion zone boundary conditions and halo profile models very important. These 
ambiguities make predictions of anti-proton and anti-deuteron spectra from dark matter 
annihilations difficult to reliably calculate. 

5 Gamma- Rays From LZP Annihilations Near the Galac- 
tic Center 

LZP annihilations can produce high-energy gamma-rays through the hadronization and de- 
cay of pions generated in the cascading of their annihilation products. In addition to this 
continuum emission, LZP annihilations can also produce 77 and 7Z final states, for example 
from box diagrams with a X s along one side and t along the other three. Such processes 
would yield the very distinctive feature of mono-energetic gamma-ray lines, although with 
a much smaller cross section than for continuum emission. The full cross sections for these 
processes have not yet been calculated for this model. For these reasons, we will focus on 
continuum gamma-rays in this section. 

Unlike positrons and other charged particles, cosmic gamma-rays travel in straight lines 
and can travel greater distances without energy loss. For these reasons, regions of the Galaxy 
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Figure 8: The cosmic positron fraction from LZP annihilations for a mass of 400 GeV 
annihilating to top quark pairs compared to the error bars projected for the PAMELA 
(larger, red) and AMS-02 (smaller, blue) experiments. We have used an annihilation cross 
section of 3 x 10~ 26 cm 3 /s, a local dark matter density of p(local)=0.4 GeV/cm 3 and a boost 
factor of 10. The solid line is the background-only prediction. 

which are expected to have higher concentrations of dark matter (and higher annihilation 
rates) can be studied with gamma-ray telescopes without integrating over the entire Galactic 
volume. In particular, the Galactic center is a promising region for such observations. 
The 7-ray flux from dark matter annihilations near the Galactic center is given by 



Here, ip is the angle between the line-of-sight (los) and the Galactic center, (av) is the LZP 
annihilation cross-section and p(r) is the dark matter density at distance, r, from the Galactic 
center. dN 1 /dE 1 is the spectrum of continuum gamma-rays produced in each annihilation. 

This expression can be further separated into two factors, one specifying the particle 
physics model (mass, cross-section and fragmentation spectrum) and another describing the 
dark matter distribution. Normalizing to the distance to the Galactic center and the local 
halo dark matter density, the latter factor can be written as: 



Defining J(AO) as the average of J{ip) over the solid angle AQ (centered on if) = 0), we 
can write 
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Figure 9: The gamma-ray spectrum from LZP annihilations in the Galactic center region. 
The dotted, dashed and solid lines are for LZP masses of 50, 400 and 800 GeV. Annihilations 
over an angular cone of size AQ ~ 5 x 10~ 5 sr were included. A Navarro-Frenk- White halo 
profile has been use as well as an annihilation cross section of 3 xl0~ 26 cm 3 /s. 

For a standard Navarro-Frenk- White (NFW) halo profile, averaged over a solid angle of 
Afl ~ 5 x 10~ 5 sr, we calculate a value of J(5 x 10~ 5 sr) ~ 5.6 x 10 3 . This can vary- 
considerably with halo model, however. For example, with a Moore et al profile [44], we 
calculate a much larger value of J(5 x 10~ 5 sr) ~ 1.9 x 10 6 . If the effects of adiabatic 
compression or accretion onto the central supermassive black hole are taken into account, 
even larger values can be found [45]. 

The spectrum of gamma-rays predicted from LZP annihilations is shown in figureEl This 
result is not substantially different than what would be expected from the annihilation of 
neutralinos or other WIMPs. Also, recent observations of a gamma-ray spectrum extending 
above ~10 TeV may pose serious challenges to identifying any potential gamma-ray signal 
from WIMP annihilations [46], assuming that this is not in fact the annihilation signature 
of a very heavy WIMP [47]. For these reasons, we do not further discuss the gamma-ray 
indirect detection channel here. 

6 Comparisons to SUSY and UED Models 

In this section, we briefly summarize the phenomenological features relevant to dark matter 
detection for an LZP Dirac RH neutrino and compare these features to those of the lightest 
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neutralino in supersymmetric models and the first Kaluza-Klein excitation of the hypercharge 
gauge boson in models with Universal Extra Dimensions. 

Mass Range: LZPs can generate the observed quantity of dark matter thermally in two 
mass ranges: near the Z-resonance (m-^zp ~ 35-50 GeV) and for considerably heavier masses 
(^lzp ^ several hundred GeV) [14, 15]. Neutralinos can span a very wide range of masses, 
although in the models most often studied, they fall in the range of 50-1000 GeV. The LKP 
in UED models requires a mass in the range of 700-900 GeV in order to generate the observed 
thermal relic density, unless other KK modes participate in the freeze-out process [3]. If this 
is the case, somewhat smaller masses are possible. A recent analysis, taking into account 
the effects of second level KK modes, indicates that the upper edge of this mass range is 
favored [4]. In any case, we should keep in mind that the precise prediction of the LKP 
relic density depends on the particular KK mass spectrum which is used and is somewhat 
model-dependent. 

Annihilation Cross Section: If any of these dark matter candidates generate the observed 
relic density thermally, they must have an annihilation cross section of (av ) ~ 3 x 10~ 26 
cm 3 /s evaluated at the freeze-out temperature. For the purpose of indirect detection, we 
are interested in the annihilation cross section in the low velocity limit, which does not 
necessarily equal the value at freeze-out, however. For most of the parameter space in this 
model, these two cross sections are not considerably different. In UED models, they are 
always the same. For neutralinos, this issue depends on the composition of the LSP and the 
spectrum of other supersymmetric particles. For the case of neutralino annihilation through 
s-channel pseudoscalar higgs exchange to fermion pairs, as if often the dominant channel 
for bino-like LSPs, the cross section at freeze-out will be very similar to the value at low 
velocities. 

Dominant Annihilation Channels: The leading annihilation channels for LZPs are shown 
in Fig. El Below 100 GeV, annihilations to light quarks, neutrinos and charged leptons are 
important. At higher masses, ti, W + W~ and hZ channels dominate. In UED models, about 
20% of LKP annihilations go to charged lepton pairs of each family. About 1.2% go to each 
flavor of neutrino pairs. The remaining annihilations go to quarks and higgs bosons. 

For neutralinos, again the composition and SUSY spectrum are important in predict- 
ing the dominant annihilation modes. For a bino-like LSP, annihilations to bb, t + t~ and, 
if kinematically allowed, ti usually dominate. For winos and higgsinos, annihilations to 
gauge and higgs bosons are often important. Since neutralinos are Majorana fermions, their 
annihilations to light fermions are chirality suppressed, unlike for the LZP and LKP. 

Elastic Scattering With Nucleons: LZP capture in the Sun can be dominated by spin- 
dependent scattering with hydrogen nuclei (protons) or by scalar scattering with helium 
nuclei. The size of these cross sections can vary a great deal, depending on the choices of 
parameters. In particular, the scale Mkk (mass of KK gauge bosons) plays a very important 
role. The bounds on the spin-independent cross section placed by the CDMS experiment [48] 
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are often exceeded for smaller values of M KK ~ 3-4 TeV [14, 15]. On the other hand, values 
as small as ~ 10~ 9 pb are possible for M KK ~10 TeV. 

For neutralinos, these cross sections can vary from in excess of the CDMS bound, to as 
small as ~ 10~ 15 pb, although in constrained models, they often fall in the range of ~ 10~ 6 - 
10~ 9 pb. In UED models, spin-dependent scattering dominates LKP capture in the Sun, 
with cross sections of 1CT 5 -1CT 6 pb. 

7 Discussion and Conclusions 

In this article, we have studied the prospects for the detection of a stable, right-handed, 
Dirac neutrino which can emerge in grand unified models in a Randall-Sundrum spacetime. 
After introducing into the model a Z 3 symmetry necessary to ensure proton stability, the 
Lightest Z 3 charged Particle (LZP) is stable and a viable dark matter candidate. We have 
studied the indirect detection prospects for such a particle through three channels: neutrinos 
produced through annihilations in the Sun, positrons produced through annihilations in the 
Galactic halo and gamma-rays produced through annihilations in the Galactic center. 

The prospects for detecting high-energy neutrinos from the Sun in this model are very 
encouraging. For a light LZP in this model (mLZP ^, 100 GeV), annihilations often produce 
neutrinos directly. The currently operating AMANDA experiment may already be sensitive 
to some of these models. In the future, the ANTARES experiment will further probe into this 
region. For heavier LZPs, annihilations to top quark pairs produce neutrinos efficiently in 
their decays. The kilometer-scale experiment IceCube, currently under construction, should 
be sensitive to most of these models in which the Kaluza-Klein scale, Mkk, is lower than 
~ 6 TeV (a favorite range as far as fine-tuning in the electroweak sector is concerned). 

Lighter LZPs also generate positrons very efficiently A 35-50 GeV LZP annihilating in 
the galactic halo produces a spectrum of positrons which matches the spectrum observed 
by the HEAT experiment quite well, and with more moderate annihilation rates than for 
neutralino dark matter. At an energy equal to the LZP mass, a discrete feature appears 
in the positron spectrum which, if identified in future experiments such as PAMELA and 
AMS-02, would provide an unambiguous signature of dark matter annihilations. 

Finally, LZP annihilations near the Galactic center may provide an observable flux of 
gamma-rays. We have calculated the spectrum of gamma-rays in this model and find that 
is it not considerably different than for the case of annihilating neutralinos. 

These results could be generalized to other Dirac RH neutrino dark matter particles, 
though it is not clear in which other contexts such a particle would naturally arise. In 4D, 
this would require the addition of vector-like neutrinos very weakly coupled to SM gauge 
bosons but interacting via an extra U(l). The coupling of SM fermions to this extra U(l) 
should be suppressed to keep their coupling to the Z (due to Z-Z' mixing) almost unshifted. 
It is also unclear what would make this particle stable in this context. Note, however, that 
via the AdS/CFT correspondence, our Randall-Sundrum scenario is dual to a 4D composite 
Higgs scenario, in which the unification of gauge couplings has recently been studied [49] . In 
this case, the LZP maps to some low-lying hadron at the composite scale. We also point out 
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that in Refs. [14,15], the strong coupling scale is close to the curvature scale so that 0(1) 
variations in our calculations are expected. These results should therefore be considered as 
representative rather than a complete description. 

Annihilations can vary from one Dirac RH neutrino dark matter model to another, de- 
pending on whether, at large v' R mass, annihilations take place via s-channel Z' exchange 
only or also via a t-channel X^-type gauge boson. On the other hand, the elastic scattering 
cross section is mainly model- independent (determined by the v' R — Z coupling). 
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